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ABSTRACT 

Aims. Employing photometric rotation periods for solar-type stars in NGC 1039 [M 34], a young, nearby open cluster, we use its mass- 
dependent rotation period distribution to derive the cluster's age in a distance independent way, i.e., the so-called gyrochronology 
method. 

Methods. We present an analysis of 55 new rotation periods, using light curves derived from differential photometry, for solar type 
stars in the open cluster NGC 1039 [M34]. We also exploit the results of a recently-completed, standardized, homogeneous BVIc 
CCD survey of the cluster, performed by the Indiana Group of the WIYN open cluster survey, in order to establish photometric cluster 
membership and assign B-V colours to each photometric variable. 

We describe a methodology for establishing the gyrochronology age for an ensemble of solar-type stars. Empirical relations between 
rotation period, photometric colour and stellar age (gyrochronology) are used to determine the age of M34. Based on its position in 
a colour-period diagram, each M 34 member is designated as being either a solid-body rotator (interface or I-star), a differentially 
rotating star (convective or C-star) or an object which is in some transitory state in between the two (gap or g-star). Fitting the period 
and photometric colour of each I-sequence star in the cluster, we derive the cluster's mean gyrochronology age. 
Results. Of the photometric variable stars in the cluster field, for which we derive a period, 47 out of 55 of them lie along the loci 
of the cluster main sequence in V/B-V and V/V-I space. We are further able to confirm kinematic membership of the cluster for half 
of the periodic variables [21/55], employing results from an on-going radial velocity survey of the cluster. For each cluster member 
identified as an I-sequence object in the colour-period diagram, we derive its individual gyrochronology age, where the mean gyro 
age of M 34 is found to be 193 ± 9 Myr. 

Conclusions. Using differential photometry, members of a young open cluster can be easily identified in a colour-magnitude diagram 
from their periodic photometric variability alone. Such periodicity can be used to establish a period-colour distribution for the cluster, 
which for M 34, we have used to derive its gyrochronology age of 193 ± 9 Myr. Formally, our gyro age of M 34 is consistent (within 
the errors) with that derived using several distance-dependent, photometric isochrone methods (250 ± 67 Myr). 

Key words. Methods: data analysis - Stars: starspots - Galaxy: open clusters and associations: NGC 1039 - Stars: fundamental 
properties (age) 



1. Introduction 



Over the past century, as instrumentation and detector technol- 
ogy have advanced, the methods by which stellar age is deter- 
mined have also evolved. While we still rely heavily upon the 
use of theoretical model isochrones, fitting empirical luminosity- 
temperature data (on an H-R diagram - e.g., Sandage 1958; 
Demarque & Larson 1964; Demarque & Gisler 1975, Mengel et 
al. 1979; and more latterly, Meynet, Mermilliod & Maeder 1993; 
D'Antona & Mazzitelli 1997; Naylor 2009), several alternative 
techniques for determining stellar age now exist which exploit 
a more diverse assembly of fundamental stellar properties, such 
as magnetic activity, elemental abundances, white dwarf cooling 
time-scales and angular momentum content. 



Inspired by the pioneering efforts of Wilson (1963), obser- 
vations of solar type stars have shown that there is a clear age 
dependence in the strength of magnetic activity indicators (eg., 
Ca II H & K and X-rays), whose emissions are presumably due 
to solar-like magnetic field activity (eg., Wilson 1964; Wilson 
& Skumanich 1964; Wilson 1966; Skumanich 1972; Barry et al. 
1981;Noyesetal. 1984;Maggio etal. 1987; Soderblom, Duncan 
& Johnson 1991; Henry et al. 1996; Mamajek & Hillenbrand 
2008). While there is not a one-to-one relationship between ac- 
tivity and age for any given star, current data samples unequiv- 
ocally show there to be a strong correlation between decreasing 
magnetic activity and increasing age. 

However, even such a well-observed phenomenon does not 
deliver unequivocal results, as illustrated for instance by the fac- 
tor of two range in ages derived by Giampapa et al. (2006) for 
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photometric and kinematic members of the 4Gyr-old open clus- 
ter M 67. Moreover, stellar activity cycles, akin to those occur- 
ring in the Sun, can also lead to difficulties in establishing ages 
for individual stars. For instance, Soderblom et al. (1991) show 
that errors on chromospheric ages are roughly 50%. Presumably, 
this error could be driven down statistically in a cluster by mea- 
suring many stars, and by removing close binaries or other con- 
taminants. Soderblom, Jones & Fischer (2001) have provided 
chromospheric Ha measurements for a significant number of 
M 34 stars, but to our knowledge, a chromospheric age for M 34 
has not yet been published. 

More recently, age determinations for young open clusters 
have been advanced based on an analysis of their light element 
abundance distributions, specifically using lithium as a tracer. 
This technique, the so-called lithium depletion boundary method 
(Basri, Marcy & Graham 1996; Rebolo et al. 1996; Stauffer, 
Schultz & Kirkpatrick 1998; Barrado y Navascues, Stauffer & 
Patten 1999), does however have both merits and limitations. 
Line strengths and elemental abundance ratios are fairly straight 
forward to measure, and are indeed fundamental stellar proper- 
ties, perhaps more so than compared to magnetic activity indica- 
tors. Unfortunately, the strength of this technique lies in detect- 
ing lithium in M-stars and cooler, which for any group of stars 
at a given distance are among the faintest. Even for nearby open 
clusters, detecting the lithium boundary thus requires use of pre- 
cious 8-10m class telescope time. Perhaps the ultimate limita- 
tion for this technique is that it is only discernibly sensitive to 
stars of $250 Myr (e.g., Stauffer et al. 1998), severely limiting 
its widespread applicability. 

The use of heavy-element abundances to establish stellar age 
has also been attempted, however its implementation is as yet 
limited. Cayrel et al. (2001) advocate using the detection of the 
singly-ionized uranium 238 line at 3859. 6A as an age proxy, al- 
though as they themselves note, the accuracy of this technique 
is currently limited by deficiencies in the nuclear data - e.g., 
poorly-known oscillator strengths. Even with the advent of well- 
calibrated atomic data for uranium species, dating of open clus- 
ters using detection of stellar uranium will only be feasible for 
clusters older than the Hyades [600Myr], due to the prohibitively 
long half-lives of the two dominant naturally occurring uranium 
species (n/2 = 7.04 x 10 8 & 4.47 x 10 9 yrs for U 235 & U 238 
respectively - CRC 2005). 

Spectroscopic and photometric observations of white dwarfs, 
in the Galactic field and in globular and open clusters, can also 
yield an age estimate of the their stellar content. The detailed 
physics of how white dwarfs cool as they age is now well de- 
veloped (Mestel 1952; Cox & Giuli 1965; Beaudet, Petrosian & 
Salpeter 1967), and their model-dependent cooling time-scales 
are well described using theoretical stellar models (Koester 
& Schonberner 1986; D'Antona & Mazzitelli 1989; Iben & 
Laughlin 1989; Wood 1992; Fontaine, Brassard & Bergeron 
2001). This technique has been historically rather difficult to 
implement in open clusters. First, even for the nearby clusters, 
their white dwarf members are the faintest objects comprising 
the mass function, typically being fainter than V=20th magni- 
tude (e.g., Rubin et al. 2008), which is challenging for all but 
the largest aperture terrestrial telescopes. Second, white dwarfs 
are intrinsically blue objects (typically B-V fl < 0.2), which until 
recently was problematic for capturing photons efficiently using 
astronomical, red sensitive spectrographs and detectors. Third, 
and most importantly, to properly derive the age of the clus- 
ter one must observe and characterize objects extending right 
down to the bottom of the white dwarf sequence; i.e., those 
objects which are the dimmest, the coolest and hence the old- 



est white dwarfs in the cluster. Consequently, in light of such 
procedural difficulties, relatively few open clusters have white 
dwarf cooling-time age estimates available in the literature (e.g., 
Richer et al. 1998; Bedin et al. 2005; Kalirai et al. 2003, 2005; 
Rubin et al. 2008). 

In terms of rotation, Kawaler (1989) uses theoretical stel- 
lar spin-down models (from Kawaler 1988) to analytically de- 
rive a relationship uniquely connecting a star's B-V photometric 
colour and age to its period. Using a photometric colour, rota- 
tion period dataset for G & K-stars in the Hyades open cluster, 
he shows the Hyades rotation age to be 4.9+ 1 . 1 x 10 8 years. More 
recently, Barnes (2003, 2007 - hereafter B03, B07 respectively) 
exploit the morphology of photometric colour-rotation period 
distributions in coeval samples of stars (e.g., open clusters and 
binary systems) to establish stellar age. In essence, this so-called 
gyrochronology technique allows one to construct rotational 
isochrones, in order to trace the boundaries of age dependent 
colour-period distributions. Crucially however, these boundary 
definitions further suggest the identification and coupling-state 
of each star's basic internal stratification in terms of their radia- 
tive zone and convection envelope. 

Solar-type stars lying close to or on the /, or interface se- 
quence, are probably rotating as solid bodies or close to it, 
and their spin-down evolution is Skumanich-like (Skumanich 
1972), i.e., directly proportional to the square root of stellar 
age (spin-down oc f 1 ^ 2 ). The most rapidly rotating solar-type 
stars at a given mass lie on or close to the C, or convective se- 
quence. These stars are thought to be rotationally de-coupled; 
that is to say, it is likely that only their outer convective envelope 
is spinning down, with an exponential time dependency (spin- 
down oc e^ w ). Early incarnations of rotational evolution mod- 
els incorporating de-coupled stellar structure were referred to as 
core-envelope (de)coupling models (c.f., Endal & Sofia 1981; 
Stauffer et al. 1984; Soderblom et al. 1993a, Jianke & Collier 
Cameron 1993). Those stars which are in a transitory state be- 
tween the C and /-sequences, representing a scenario whereby 
the outer convective envelope is re-coupling, probably magneto- 
hydrodynamically, to the inner radiative zone, constitute the so- 
called gap or g-stars in the gyrochronology paradigm. 

Establishing stellar age for open clusters using gyrochronol- 
ogy models of course has its limitations. For instance for open 
clusters, enough stars must be photometrically monitored in or- 
der to derive rotation periods, so that a clear, well-defined dis- 
tribution of interface, gap, or convective sequence objects is ap- 
parent. Such an observing programme typically requires exten- 
sive allocations of telescope time and considerable human effort. 
Moreover, differential rotation of solar-type stars and multiple 
star-spot groups on their surfaces can also introduce ambiguities 
into gyrochronology analyzes, since it can act to smear-out the 
distribution of rotation periods of a given mass. 

For small samples of rotation periods in open clusters, the 
gyrochronology method is not applicable in multiple systems 
where other effects might interfere with, or even overwhelm, the 
regular wind-related loss of angular momentum. For instance, 
close binaries (a $0. 1 AU) experience a tidal torque which acts 
to synchronize their rotation with the orbital motion of the sys- 
tem, in addition to any magnetic torque they experience due to 
stellar winds (Zahn 1989; 1994); thus, as they evolve along the 
main-sequence, components in the closest binaries tend to ro- 
tate faster (on average) than single stars of similar mass and 
age. Observations of RS CVn and BY Dra systems have shown 
Zahn's theoretical framework to be correct (Hall 1976; Bopp & 
Fekel 1977; Fekel, Moffet & Henry 1986). Binary/multiple sys- 
tems must therefore be avoided, where the data allow, in estab- 
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lishing gyrochronology ages of stars (see Meibom, Mathieu & 
Stassun 2006 for a more thorough discussion). 

A recent International Astronomical Union symposium 
(IAUS 258), dedicated to the subject of The Age of Stars, now 
provides our community with a thorough overview of an histor- 
ical and modern approach to understanding the how and why 
of stellar ages. A careful perusal of the symposium proceedings 
allows, for both the novice and expert alike, a detailed insight 
into the theoretical framework unpinning stellar age determina- 
tions as well as the observational data calibrating and constrain- 
ing such models. Of special interest to this manuscript were the 
oral presentations by Barnes (2009), Jeffries (2009) and Meibom 
(2009). 

The genesis of this gyrochronology project for M 34, a =t 
200 Myr open cluster, occurred nearly ten years ago. At that 
epoch, SAB and GWL obtained a differential photometry dataset 
at Lowell Observatory over 17-nights for the Western half of 
the cluster, with the goal of deriving rotation periods for photo- 
metric variables in the field. The project lay dormant for some 
years until we noticed that Irwin et al. (2006) published new ro- 
tation period data for a sample of solar-type stars in the cluster. 
Unfortunately in terms of rotation, there are several fundamental 
problems with the Irwin et al. dataset which hamper our ability to 
pursue a rigorous gyrochronology analysis using their data (see 
§ l4.4l for details). A natural progression and indeed amelioration 
of their period distribution analysis is to derive the gyrochronol- 
ogy properties of this cluster, now incorporating rotation periods 
derived from the differential photometry that we had obtained 
during the 1998 Lowell campaign. In this manuscript therefore, 
we now report an analysis of our differential photometry dataset 
for the cluster. 

New BVI photometry is reported in §[3]for each photometric 
variable discovered in our Lowell observing campaign. Rotation 
periods derived from light-curve analysis of the photometric 
variables in our field of view are described in §|4] which includes 
a comparison with the Irwin et al. study. Using a clean sample 
of rotation periods for photometric and/or kinematic members of 
the M 34 cluster, we describe in §[5]how colour-period data are 
used in undertaking a gyrochronology analysis, which results in 
the derivation of the rotation age of the cluster. Also included 
is a commentary on how error in the gyrochronology age can 
be ascribed to disparate physical properties of the colour-period 
dataset, such as cluster non-membership, binarity and differen- 
tial rotation. 



2. The Open Cluster Messier 34 

M 34 (NGC 1039) is a young (^200 Myr; see TableQ) open clus- 
ter on the border of Perseus and Andromeda [RA(2000): 02 A 42 m , 
DEC(2000): +42°46']. It is a 470 pc away, out of the Galactic 
plane [b a —16°], and has a reddening of Eb_i/=0.07 (Canterna, 
Perry & Crawford 1979). These characteristics suggested to us 
that it would make a good target for rotation period determina- 
tions, including a derivation of its rotational age via the method- 
ology of gyrochronology. 

2.1. Extant Observations ofM34 

Prior work on M 34 includes photoelectric UBV observations for 
57 stars bluer than F3 in the M 34 field by Johnson (1954) and 
of 43 similar stars by Cester et al. (1977). Canterna et al. (1979) 
reported photoelectric photometry in Stromgren uvby & HB, for 
some 42 M34 stars. Jones & Prosser (1996) published a more 



extensive and deeper (down to V=16.2) survey using BVI CCD 
observations. An even deeper (V^ 24) VI CCD survey has been 
recently published by Irwin et al. (2006). Both of the latter two 
studies have certain irregularities with the standardization proce- 
dures that are evident in the papers. Consequently, we use stan- 
dardized photometry from a new BVRI survey by the Indiana 
group of the WIYN open cluster study (Deliyannis et al. in prep) 
in our analysis (see §|3j. 

Two proper motion surveys have been conducted in M 34. 
The first, by Ianna & Schlemmer (1993), has a limiting mag- 
nitude of V=14.5, and the second, by Jones & Prosser (1996), 
one of V=16.2. Reliable membership probabilities from these 
studies are limited to stars brighter than V=13.5 and V=14.5, 
respectively, and although >80% of the stars in our study are 
fainter than either of the proper motion surveys, we use them 
where available. Jones et al. (1997) have provided radial veloc- 
ities for some 47 cluster members but few of these overlap with 
our photometric period sample. However, we take recourse to an 
on-going multi-epoch WIYN HYDRA radial velocity survey of 
the cluster (Meibom et al. in prep) to provide a reliable kinematic 
membership criterion for nearly half (25/55) of the M 34 stars for 
which we report new photometric periods. A thorough discus- 
sion of extant membership probabilities for M 34 stars compris- 
ing our gyrochronology sample, as well as new high resolution 
spectroscopic results for two photometric variables in our sam- 
ple, is presented in Appendix lAl 

Other meritorious work on M 34 includes a study of the clus- 
ter white dwarfs by Rubin et al. (2008), who derive a minimum 
age of the cluster to be 64 Myr. Although Rubin et al. were not 
explicitly attempting an age derivation, such a result must be 
considered with some caution because of the lack of surety of 
cluster membership for their targets, the small sample size (5 
objects), and the fact that only the very top of the cluster white 
dwarf sequence was sampled, whereby the coolest, and hence 
oldest white dwarfs in the cluster were not observed. 

Projected equatorial rotation velocities [v sin;] for the 
earlier-type, higher mass stars in M 34 are provided by Ianna 
(1970), Hartoog (1977), and Landstreet et al. (2008), whereas 
for its solar-type members velocities are published by Jones et 
al. (1997) and Soderblom et al. (2001). For stars in common, 
we provide a comparison between these v sin i data and our new 
photometric periods in Appendix [B] By removing the dual am- 
biguities in v sin i data (the obvious sin i factor, and the stellar 
radius), photometric rotation period studies far and above super- 
sede v sin i studies - provided of course that the rotation periods 
are measured correctly. 

Irwin et al. (2006) have provided rotation periods for some 
105 stars in M34, using photometric data obtained under con- 
siderably more restrictive conditions (observing window, and fil- 
ters) than ours. Although half the stars in common between their 
study and our Lowell sample yield essentially the same period 
from multiple time-series, some discrepancies exist. A detailed 
comparison between their periods and ours is presented in § 14.41 

3. New Indiana CCD photometry 

A new, comprehensive BVRI survey of star fields in the direc- 
tion of the M 34 cluster is now available from the Indiana Group 
of the WIYN Open Cluster Survey (Deliyannis et al. in prep). 
Photometric observations of the central 40 x 40 of the clus- 
ter, including a 100 Landolt (1992) standard stars, were ob- 
tained using the WIYN 0.9-meter telescope during the night 
of 24-Oct-1998. The 2048 x 2048 pixel T2KA Tekronics CCD 
was employed at the f/7.5 Cassegrain focus, which with its 
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'.'68 pixel 1 plate scale yields an effective field of view of 
s 23 x23 . Multiple exposures in each filter were obtained over a 
2x2 mosaic for the cluster. Images were analyzed using the dedi- 
cated DAOPHOT II photometry data reduction package. Landolt 
standard star frames were processed using aperture photome- 
try, whereas object frames were fitted with spatially-variable 
point spread functions [PSFs], and spatially-variable aperture 
corrections. Both the PSFs and aperture corrections were de- 
termined empirically using hundreds of isolated stars on each 
frame. External errors in placing instrumental magnitudes onto 
the standard system are < 0.02 magnitudes. 

The results of the Indiana Group's photometric survey for 
stars relevant to this study are presented in Table|2] as well as co- 
ordinates for each star derived from SuperCosmos survey images 
(Hambly et al. 2001a, b; Hambly, Irwin & MacGillivray 2001). 
We also provide a kinematic membership assessment from an 
on-going survey of the cluster (Meibom et al. in prep). For two 
M 34 targets for which we derive a period, F4_1925 & F4_2226, 
there are no fi-band photometric data available. In this instance, 
we rely upon a field star V-I to B-V relationship (Caldwell et 
al. 1993), and an assumed E B _y=0.07 (Canterna et al. 1979), to 
transform these stars' V-Ic colours into B-V ones. 



4. New Photometric Periods in M34 

4.1. Observational Programme 

We initiated a differential photometry programme of observa- 
tions for solar-type stars in M34 during September-October 
1998, consisting of four separate sequences of CCD images in 
the V and I filters, acquired respectively during a 17 night ob- 
serving run (UT980929 through 981015) at Lowell Observatory, 
Flagstaff, Arizona, USA. The dataset consists of ~ 60 visitations 
to each of two 19'xl9' fields of the Western region of M34. 
Each visitation typically consisted of short (60, 40 sec) and long 
(600, 400 sec) exposures in the V and I filters respectively, us- 
ing the 2048x2048 SITe chip mounted at the Cassegrain focus 
of the 42inch Hall telescope of Lowell Observatory. Each field 
received 4-7 visitations per night in each filter, over the first 12 
nights, barring UT981008 which was lost to clouds, and one vis- 
itation each on the remaining 5 nights. The time-series analysis 
is performed on the instrumental system with respect to chosen 
reference images on UT981001. 

4.2. Data Analysis and Time Series Analysis: 

The individual images were de-biased and flat-fielded using 
standard procedures in the IRAF3 image reduction utility. 
Photometry for all objects on the processed images was com- 
puted via PSF fitting using DAOPHOT II and ALLSTAR II 
(Stetson, Davis & Crabtree 1991). Typically, we have used 50 
PSF stars per image and a quadratically-varying PSF. The pho- 
tometry from the individual frames was cross-matched and con- 
verted into time series using the DAOMATCH/DAOMASTER 
routines from Stetson (1992). These routines match stars ap- 
pearing in multiple images according to user-supplied specifi- 
cations, calculate the photometric offsets between frames, and 
output the corrected time series. The uncertainty in the frame-to- 
frame magnitude offsets was less than 0.001 mag for the frames 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



that were eventually used because of the large number (^1000) 
of stars matched from frame to frame. This procedure left us 
with totals of ~ 4000, ~ 3000, ~ 2700, and ~ 2200 objects re- 
spectively with time series in each of the long-I, long-V, short-I 
and short-V exposures for the western M 34 fields. 

We tested every object in the field for variability without re- 
gard to its position in a colour-magnitude diagram, or any other 
consideration except its intrinsic variability. Also, we performed 
the variability/periodicity analysis independently in each expo- 
sure time/filter combination, in order to confirm or deny it, and 
to calculate period error. The following procedure was adopted 
for each time-series file. The time series files for each field were 
sigma-clipped to eliminate outlierfl and then we calculated a 
reduced chi square (x 2 ) for each time-series, to glean out the 
most variable objects. Every object with x 1 > 3 was subjected 
to a periodicity search using the CLEAN algorithm of Roberts, 
Lehar & Dreher (1987). The results were examined visually, 
and phased about reliable peaks in the power spectrum to yield 
phased light curves for the periodic variables. Finally, results of 
the (up to) four independent time-series analyzes on each object 
were merged together to confirm or deny the periodicity. For 
these stars, we calculate period errors based on multiple mea- 
surements. 

As a result of the foregoing analysis, we have been able to 
identify forty-seven (47) objects with periods confirmed in at 
least two filter, exposure combinations. In addition, we derive 
periods for eight (8) objects for which we can only detect pe- 
riodicity in one filter/exposure. The period analysis results are 
listed in Table|3] which also contains details of their derived peri- 
ods, and standard deviations of magnitude variation in each time- 
series. Stars having only one period determination are assigned 
period errors based on the widths of the peaks in the power spec- 
tra. Phased photometric light-curves, on the instrumental system, 
are presented for each of our M 34 variables in Appendix ICl 

Our 17-night observing window allows us to detect two or 
more phases for all periods up to about 8.5 days. However this 
cut is probably too conservative, and because the slow rotators 
are of great interest to us, we extend our analysis to periods 
longer than 8.5 days. We are confident about periods as long 
as three quarters of the 17d baseline, i.e. ~ 13d; however, the pe- 
riod errors are larger for such stars. For instance, F3J3469 (field 
3) has a period of ~ 12. 9d in all 4 exposure time/filter time- 
series. Periods longer than our 1 3-day comfort timescale abut the 
sensitivity limit of our observing window, which result in conse- 
quentially larger period errors for the [2/55] systems for which 
we report periods >13 days (see Table[3]). 

4.3. Establishing Membership of M34 from Photometric 
Variability Alone 

Constructing a colour-magnitude diagram is the simplest way 
of establishing possible membership of an open cluster because 
stars located far from the cluster sequence are guaranteed to 
be non-members. Stars on the cluster sequence may or may 
not be bona fide cluster members, depending on the extent of 
fore- and background contamination from the Galactic field. For 
stars forming an apparent cluster main sequence, additional in- 
formation is needed to distinguish between true members and 
field-star, non-member interlopers. Classical, well established 



2 This procedure lowers our sensitivity to eclipsing binaries and flares 
or other one-time variables, but has the virtue of eliminating occasional 
bad photometric points, which could reduce the efficacy of the period- 
searching software for rotational variables. 
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methods for member, non-member triage involve astrometric 
(proper motion) measurements or measuring radial velocities 
(e.g., Stauffer 1984, Meynet et al. 1993, Perryman et al. 1998, 
Raboud & Mermilliod 1998, Mermilliod et al. 2008), however 
such measurements are not always practical to acquire. In their 
absence, magnetic activity (e.g., X-rays) or lithium abundance 
indicators can be used to distinguish between cluster and non- 
cluster stars (e.g., Prosser et al. 1995a; James & Jeffries 1997; 
Cargile, James & Platais 2009). Here, we propose that photo- 
metric periodicity can also be exploited to easily establish open 
cluster membership in a manner unbiased with respect to bright- 
ness, mass (at least for solar-type stars) and distance. 

In Fig.,[T] we plot colour magnitude diagrams of the cluster 
in BV and VI colours, along with D' Antona & Mazzitelli (1997) 
theoretical isochrones for 250 Myr & 470 pc in both colours. The 
55 periodic variables in our Lowell fields are also identified on 
the figure. Photometric and/or radial velocity members, 43 in 
number, are flagged with red squares, whereas non-members, 12 
in number, are flagged with blue circles (their periodicity and 
other information are tabulated in Tables [2] & [3). The cluster 
sequence is traced out very well by the periodic variables. In 
fact, in this particular case, photometric periodicity information 
alone provides a 78% (43/55) probability of cluster member- 
ship. We believe the reason for this behavior is that the cluster 
stars are younger and hence far more active and variable than 
the mostly older disk stars in the field. One might even envi- 
sion using photometric variability to trace out cluster sequences 
in sparsely populated young clusters or background-dominated 
cluster fields. 

4.4. Comparison with prior periods 

Irwin et al. (2006) have published a prior study of periodic vari- 
ables in M 34, and before proceeding further, we compare our 
periods with theirs for stars in common between the two stud- 
ies. This comparison is graphically represented in Fig. ,[2] where 
symbols are assigned to stars having a varying number of inde- 
pendent period derivations (1-4) in our dataset. Of the 19 stars in 
common between the Irwin et al. sample and our own Lowell 
campaign, roughly half (9/19) have periods that agree within 
10%. For all but two of the remaining 10 stars, Irwin et al. peri- 
ods are significantly shorter than ours, where the origins for such 
discrepancies are not entirely clear to us. 

Through comparison with an as yet unpublished photomet- 
ric variability study of M 34 we can establish confidence in the 
Lowell campaign periods, and advocate their preferential use 
compared to the Irwin et al. ones. Three stars in common to the 
Irwin et al. study and our Lowell campaign, F3J3413, F4J3234, 
and F4_l 147, have Meibom et al. (in prep) periods, whose pho- 
tometric variability was characterized over several months worth 
of observation data (e.g., see Meibom et al. 2009 for details 
of a similar observing programme for the M 35 cluster). These 
data, flagged stars in Fig.,|2]and listed in Table|4] show that the 
Meibom et al. periods are an excellent match to the Lowell pe- 
riods (to within 1%), whereas the Irwin et al. periods are poorly 
matched. 

Proving the correctness of one period derivation over another 
is not an easy task. However, we note that our dataset has a 17- 
night baseline compared with the 10-night baseline of Irwin et 
al., suggesting that we have better sensitivity to longer period 
systems. Our data were also acquired all night, rather than re- 
stricted to half nights as the Irwin et al. observational data were. 
These facts strongly suggest that our data are less sensitive to 
period aliasing; indeed, several of the Irwin et al. stars lie on or 



close to particular alias curves, as shown in Fig., [2] Finally, we 
note that most of our periods [17/19] are confirmed in at least 
a second filter/exposure time combination in the Lowell cam- 
paign data. Consequently, in light of the concerns that we have 
for aliasing effects in some of the Irwin et al. data, we will pref- 
erentially use our own periods for the remainder of this work. 

5. The Photometric Period Distribution in M34 

Our rotation period distribution for photometric and/or kine- 
matic members of M 34, as a function of intrinsic B - V and 
V - I colours, is displayed in Fig., [3] In agreement with prior 
results in various open clusters (B03, B07), there are no stars lo- 
cated in the upper left portion of the colour-period diagrams. The 
Sun-like [(B-V)„ - 0.64] stars in our M 34 sample are rotating 
about an order of magnitude faster than the Sun. However the 
lower mass stars, mid-K and later, have an order-of-magnitude 
spread in rotation period, as has been seen before in other young 
open clusters. Further detailed interpretation of the M 34 colour- 
period diagram requires a comparison of the colour-period dis- 
tributions for other young clusters. 

Rotation period distributions for two comparison clusters 
whose ages straddle that of M 34 will allow us to critically in- 
vestigate the gyrochronological properties of M 34. To this end, 
we display colour-period diagrams for the ^135 Myr M35 and 
600 Myr Hyades clusters in Fig., [4] Colour-period data are taken 
from the recent study of M 35 by Meibom et al. (2009), and for 
the Hyades, from Radick et al. (1987) and Prosser et al. (1995b). 
This side-by-side comparison shows that the ~135 Myr-old M 35 
cluster has two distinct sequences of rotating stars indicated 
by the solid and dashed line in the figure, while the older 
Hyades cluster has only one. B03 called these sequences T (for 
Interface) and 'C (for Convective), with 'g' (for gap) stars in 
transition from C to I (see § Q] for a discussion of the angular 
momentum characteristics of Cgl stars). What we see plainly 
from the comparison is that by the Hyades age, most rotating 
solar- type stars have been transformed into I-type stars. Graphs 
of this C-g-I transformation for a series of clusters can be found 
in B03 and Meibom et al. (2009). 

We now take these gyrochronology sequences from the M 35 
colour-period distribution, and overplot them on the M 34 data 
in the left hand panel of Fig., [5] This allows us to classify the 
gyrochronology status of the (far fewer) M 34 stars roughly, as 
indicated by the Cgl symbols in the figure. We note that its 
C sequence is in the same approximate position as that of M 35, 
but its I sequence lies above that of M35, as expected from its 
older age. The manner of this inter-comparison suggests that 
only the stars shown in the right panel of the figure ought to 
be identified as I-sequence stars in M 34, which we henceforth 
do. 

5.1. A Gyrochronology Age for M 34 

It is desirable to derive an age for M 34 using our photometric 
period dataset, and the method of gyrochronology (B03, B07) 
will allow us to pursue such an analysis. This exercise will allow 
us to examine whether the gyro age derived is physically reason- 
able, and hence whether it might be usable for other clusters in 
future. The basis for such a gyro age calculation relies upon; 

1 . Deriving the true rotation period distribution of the cluster, 

2. Removing possible cluster non-members, and 

3. Choosing, and using only, the I-sequence stars in the cluster 
for the analysis. 
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For M 34, we have performed these three steps to the best of 
our ability, as described above, with the resulting set of candi- 
date I-sequence stars in the cluster being displayed in the right- 
hand panel of Fig., We note that the forthcoming study by 
Meibom et al. (in prep), which contains many more M 34 stars, 
will make the Cgl distinction in this cluster more obvious, and 
verify our classifications made here. Based on fitting the open 
cluster colour-period distributions then available, B07 suggested 
that the positions of the I sequences of open clusters, or indeed 
that of field stars, follows the relationship; 

Prof = f(B - V).g(t) (1) 

where P, or , (B-V), and t represent the measured rotation pe- 
riod, measured de-reddened colour, and the age respectively, and 
/ and g are fitted functions of the stellar colour and age respec- 
tively, where he proposed that 

f(B-V)=a(B-V-c) h (2) 
and 

g(t) = f (3) 

for which Barnes (2007) derived coefficients of a — 0.7725 ± 
0.0110, b = 0.601 ± 0.024, c = 0.4, and n = 0.5189 + 0.0070. 
However, based on the most extensive open cluster rotation pe- 
riod study to date, that of M 35 by Meibom et al. (2009 - dis- 
played in the left-hand panel of Fig.,|4]l, the coefficients of the I- 
sequence fit have been recently updated to be a — 0.770 + 0.014, 
b = 0.553 ± 0.052, and c = 0.472 ± 0.027 with a fixed n = 0.52. 
However, the index n is set by the solar calibration, by de- 
manding that the gyrochronology relationship for the I-sequence 
yields the solar rotation period at the solar age, the rationale for 
which is explained at length in B07. We note that, using their 
extensive rotation period dataset for M 35, Meibom et al. re-fit 
the colour function coefficients of Equ., 2 but did not re-fit the 
age exponent, «, of Equ., 3. Therefore to be mathematically rig- 
orous, we employ the Meibom et al. (2009) a,b & c coefficients 
and the solar data summarized in Barnes (2007) to re-derive the 
g(t) exponent to be n = 0.5344 ± 0.0015. 

Taking the Meibom et al. (2009) a,b & c coefficients and the 
newly derived n exponent, we fit the Equ., 1 relationship to the 
colour-period data of M 34 I-sequence stars, those identified in 
Fig. ,[5] in order to determine the cluster's gyro age, the one free 
parameter in the fit. The results of this procedure are reported 
in Table [5] An unweighted fit to these data yields a rotation age 
for M34 to be 193 ± 9 MyiQ (+ possible systematic errors in 
gyrochronology), where the quoted error is the error on the mean 
for these data. Such a gyro age for the cluster is consistent with 
the Jones & Prosser and Meynet et al. isochronal ages for M 34. 

5.2. Variance of the Derived Periods about the Fit 

Let us see if we can clarify the origin of items in the error bud- 
get contributing to the uncertainty in the cluster's gyrochronol- 
ogy age, which are in addition to any systematic error of the gy- 
rochronology method itself. We can estimate the magnitude of 
the error budget by considering the additive contribution of each 



3 Using all of the Meibom et al. (2009) coefficients, including n = 
0.52, a resultant gyrochronology age for M 34 of 223 ±11 Myr would 
be obtained. 



star's period variance^ We might expect contributions to the pe- 
riod variance from cluster non-members, tidally-interacting bi- 
naries within the cluster, a possible age spread within the cluster, 
period errors, differential rotation, and finally, one from initial 
variations in the cluster's natal angular momentum distribution. 
The data listed in Table [5] show that the total period variance 
about the fit is ~ 43.7 days 2 . 

It is beyond the scope of this paper to consider possible 
systematic errors in the gyrochronology method, which is cali- 
brated using both the Sun and a selection of young open clusters 
with measured photometric periods. An understanding of this 
error will probably require another decade of rotational work on 
open clusters. It would be unsurprising if such a study eventu- 
ally identified a 15-20% systematic error in the application of 
gyrochronology to real colour-period samples. 

The cluster colour-magnitude diagrams (see Fig. [TJ suggest 
that 3 periodic field variables could have sneaked into our 43 
member star sample had we not had the radial velocity measure- 
ments to reject them. This suggests that of our 22 purely photo- 
metrically chosen members (see Table [3}, one possibly two are 
non-members, making this unlikely to be any significant contrib- 
utor to the error. The mean variance contributor of any given star 
in our M 34 period sample is 1 .68 days 2 (see Table |5j, which 
means that non-member variance contribution from photomet- 
ric members is at most ^ 3.4 days 2 (~ 8% of the variance to- 
tal). Of the 26 candidate I-sequence stars in M34, twelve of 
them are classified as photometric-only members of the clus- 
ter, i.e., having no kinematic classification (with a mean period 
variance of 2.33 days 2 ). Of these twelve, only three have notice- 
ably large period variances (F4_1357, 7.84 days 2 ; F4.1007, 5.76 
days 2 ; F3J3664, 4.41 days 2 ). In any case contribution to period 
variance, and hence gyro age errors, by cluster non-members is 
at a low level. 

There is also a possibility that some scatter in the gyro ages 
for individual M34 stars is contributed by systems in tidally 
locked binaries, which will be addressed in greater detail in the 
Meibom et al. (in prep) study. We note that of the candidate I- 
sequence stars reported in Table [5] 14/24 have kinematic clas- 
sifications. Only one of these stars is a binary system, with no 
data being available as to whether it is tidally-locked or not. On 
average, one such data point would contribute 1.68 days 2 to the 
total variance, or ~ 4%. 

Observational studies of other young open clusters and as- 
sociations suggest that they are coeval to within about a few 
Myr (e.g., Wichmann et al. 2000; Palla & Randich 2005; James 
et al. 2006). While spreads in observed quantities that also 
change with age are generally believed to originate in other 
ways, we note that recent theoretical work by Baraffe, Chabrier 
& Gallardo (2009) shows that as a phenomenon, episodic accre- 
tion onto young solar-type stars can also act to introduce a lumi- 
nosity spread, and hence age spread, in the main sequence loci of 
Hertzsprung Russell diagrams. An intrinsic age spread of lOMyr 
(e.g., Jeffries 2009) at B-Vo colours of 0.60, 1.00 & 1.40 in our 
193Myr gyrochronology fit yields a maximum period variance 
of only 0.1 days 2 . Such an intrinsic internal age spread within 
the cluster will therefore not contribute any significant variance. 



4 We define each star's period variance as (Measured Period - Fitted 
Period) 2 , where the measured period is the photometric variability pe- 
riod determined from our differential photometry and the fitted period 
is that derived from Equ 1, using dereddened B-V colours, an assumed 
mean gyro age of 193 Myr and the Meibom et al. (2009) a, b, c coeffi- 
cients, with an n exponent, n = 0.5344 ± 0.0015. 
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Inspection of Fig. [5] and Table [3] shows that photometric pe- 
riod errors are quite large for the longer periods in our sample. 
We attribute this effect to the 17-night baseline of observations 
and expect that these errors will decline with a longer baseline 
dataset, such as the forthcoming Meibom et al. (in prep) sam- 
ple. Indeed, an informal comparison shows that they will be able 
to define the cluster's I-sequence better than our dataset allows. 
Another indication of this is that M34 I-sequence stars with 
smaller periods errors are, in general, located closer to the fitted 
sequence (see Fig. [§}. Based solely on the period measurement 
errors themselves (see Table [5]), the variance contribution due 
to imprecise period measurement is 15.5 days 2 (35.5% of the 
total). 

Differential rotation must be present in this sample, in the 
sense that two stars with the same equatorial velocity might 
present spot groups at observed at different latitudes, introduc- 
ing some period scatter in the observations. We do not however 
expect this to contribute much to the observed variance, based 
on the parametrization of its contribution to the fractional pe- 
riod variance via 6P/P = io _1 - 85 P°- 3 , following the analysis in 
Donahue, Saar & Baliunas (1996) and Barnes (2007). For the 
I-sequence stars listed in Table [5] this effect contributes only 
something like 2.06 days 2 to the variance, equivalent to 5% 
of the total variance budget. 

This leaves a remainder variance of 21.1 days 2 (48% of the 
total) to be accounted for by the effects of either; (a) the sys- 
tematic difference between the gyro fit and the observation data, 
or (b) the residual influence of initial variations. Because these 
two contributors have differing dependencies on cluster and ob- 
servational parameters, it should be possible to understand their 
separate contributions as the rotation period database of open 
cluster stars expands. 



6. Discussion and Summary 

We have presented the results of a 17-night differential photom- 
etry campaign over the Western half of the central region of the 
open cluster M 34. For all photometrically variable stars, we con- 
struct differential photometry light-curves, from which we derive 
periodicity. It is assumed that the photometric variability of ob- 
jects consistent with cluster membership is due to the presence of 
magnetic-field induced starspots on their surfaces, rotating with 
the star; the derived periods are therefore representative of the 
angular velocity of the star. 

In order to assess cluster membership, we exploit an exten- 
sive standardized CCD survey of the cluster, which shows that 
the majority of photometric variable stars [47/55] lie along the 
loci of the cluster main sequence in V/B-V and V/V-I space. 
Moreover, we are able to confirm kinematic membership of the 
cluster for 21 stars from an on-going radial velocity survey of 
the cluster (Meibom et al. in prep), 5 of which show radial ve- 
locity variations (multiple systems). Four (4) of the photomet- 
ric variables are kinematic non-members. Stars which are either 
photometric or kinematic non-members were excluded from the 
gyrochronology analysis. 

We note the existence of another photometric period dataset 
for solar-type stars in the M 34 cluster (Irwin et al. 2006). We 
present an analysis comparing periods for stars in common be- 
tween our Lowell campaign and their study. We detail several 
concerns that we have with the integrity of their dataset, notably 
their short observing window (< 10 nights) and their half-night 
observing cadence. After having verified that several Irwin et 
al. periods are most likely measurement aliases, and that our 



Lowell campaign periods agree those derived from an indepen- 
dent, separate, long time base-line differential photometry sur- 
vey, we preferentially employ our period results in pursuing a 
gyrochronology age analysis for this cluster. 

Rotation periods for bona fide cluster members, in concert 
with photometric B-V colours, are used to create colour-period 
distributions, from which we outline a specific methodology for 
performing a gyrochronology analysis. In order to calculate the 
actual gyrochronology age of the cluster, we assign each M 34 
member its Cgl status, as judged from its position in the colour- 
period diagram. I-sequence stars in M 34 are classified specifi- 
cally as those objects in the colour-period diagram lying above 
the I-sequence locus of the younger M35 cluster, yielding a 
mean gyro age for all I-sequence stars in M 34 of 193 ±9 Myr. 

There are two existing sets of age determination for the 
M 34 cluster, using the traditional isochrone fitting method and 
the white dwarf cooling time-scale. Together with gyrochronol- 
ogy, all three methods suffer from various forms of dependency 
on stellar models under-pinning their theoretical frameworks. 
Furthermore, both the isochrone fitting and white dwarf meth- 
ods suffer from a dependence on (or an assumption of) cluster 
distance, while gyrochronology is free from distance as an input 
parameter. 

The isochronal age of M 34 is not that well defined (see §|2]& 
Table[T]i, with a broad range in determined values. Unfortunately, 
the lack of observed giants in the cluster, and its broad turn-off 
locus in the V/B-V colour magnitude diagram precludes a pre- 
cise assessment of its isochronal age. Based on existing studies, 
the mean isochronal age of the cluster is 250 ± 67 Myr. 

The white dwarf cooling time-scale age for M34 is also 
problematic because its result is based upon only five (5) ob- 
jects whose membership of the cluster is not yet confirmed. 
Furthermore, the terminus of the white dwarf sequence has not 
yet been observed. In simple terms, the faintest, coolest, and 
therefore oldest white dwarfs have yet to be characterized, and 
the cluster's cooling time-scale age of 64.0 + 12.9 Myr (Rubin et 
al. 2008) must be considered as a lower limit. 

We conclude by noting that a study of M 34 by Meibom et 
al., including both a multi-month differential photometry cam- 
paign and also a multi-year radial velocity membership and bi- 
narity campaign, similar to their prior work in M 35, is in prepa- 
ration. Their survey will yield several hundred new photometric 
periods of M 34 solar-type stars, the I- and C-sequences of which 
should be far more clearly defined than they are from our Lowell 
campaign results. Their dataset should allow them to better de- 
fine the gyrochronology age of M 34, and we await their results 
eagerly. 
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Table 2. Astrometry, photometry and membership assignments are presented for those M 34 stars for which photometric periods are derived in the Lowell campaign. 



Internal" 


RA" 


DEC* 


JP 96 e 


yd 


Verr d B-V d B-Verr d 


V-I" V-Ierr d Phot 6 


RV e 


Identifier 


(J2000) 


(J2000) 








Mem 


Mem 



F3.0172 
F3.0176 
F3.0215 
F3.0258 
F3.0306 
F3.0320 
F3.0383 
F3.0430 
F3.0469 
F3.0485 
F3.0487 
F3.0600 
F3.0664 
F4.0136 
F4.0169 
F4.0194 
F4.0234 
F4.0303 
F4.0317 
F4.0327 
F4.0335 
F4.0450 
F4.0667 
F4.0695 
F4.0730 
F4.0736 
F4.0803 
F4.0804 
F4.0852 
F4.0942 
F4.1007 
F4.1020 
F4.1055 
F4.1147 
F4.1315 
F4.1357 
F4.1369 
F4.1392 
F4.1566 
F4.1617 
F4.1793 
F4.1925 / 
F4_2226 / 



02 41 05.141 
02 41 36.611 
02 41 00.981 
02 42 02.552 
02 41 57.472 
02 41 37.557 
02 40 49.882 
02 41 50.490 
02 40 49.090 
02 41 06.638 
02 40 44.889 
02 41 18.447 
02 41 20.454 
02 41 47.384 
02 41 04.952 
02 41 28.142 
02 41 33.429 
02 41 49.932 
02 4149.354 
02 4144.167 
02 41 35.247 
02 41 23.098 
02 41 06.187 
02 4142.922 
02 40 49.647 
02 40 48.542 
02 41 50.301 
02 4148.974 
02 41 53.246 
02 4146.371 
02 40 50.752 
02 40 36.693 
02 40 42.791 
02 4143.854 
02 40 53.997 
02 4144.198 
02 41 51.600 
02 41 38.851 
0241 31.131 
02 41 26.292 
02 41 20.042 
02 41 21.947 
02 40 31.108 



+42 56 43.202 
+42 50 20.288 
+42 52 46.695 
+42 51 51.721 
+43 00 26.530 
+42 57 21.593 
+42 55 31.520 
+42 58 10.016 
+42 48 20.648 
+42 48 22.011 
+42 53 47.221 
+42 58 21.342 
+42 58 52.278 
+42 43 38.525 
+42 46 51.343 
+42 38 37.095 
+42 42 11.625 
+42 37 13.906 
+42 36 37.613 
+42 46 07.333 
+42 41 02.323 
+42 40 15.954 
+42 46 20.783 
+42 33 13.910 
+42 46 54.799 
+42 39 25.716 
+42 44 37.371 
+42 39 59.680 
+42 35 26.394 
+42 32 32.128 
+42 41 19.844 
+42 41 32.500 
+42 38 59.153 
+42 45 07.859 
+42 38 48.513 
+42 35 36.074 
+42 29 50.864 
+42 43 01.148 
+42 41 14.283 
+42 30 15.124 
+42 39 23.652 
+42 32 32.637 
+42 40 25.105 



JP 49 
JP 167 
JP41 
JP 289 
JP 265 
JP 172 



JP213 
JP 50 
JP 133 
JP 148 
JP 227 
JP 224 
JP 199 
JP 158 
JP 113 
JP52 
JP 197 
JP 18 

JP 229 



14.584 
14.659 
14.918 
15.236 
15.503 
15.664 
16.179 
16.424 
16.604 
16.732 
16.818 
17.394 
17.862 
13.151 
13.510 
13.662 
13.959 
14.366 
14.410 
14.483 
14.546 
14.960 
15.623 
15.643 
15.741 
15.844 
16.025 
16.034 
16.294 
16.341 
16.839 
16.608 
16.639 
17.009 
17.347 
17.476 
17.496 
17.654 
18.132 
18.165 
18.606 
19.040 
19.659 



0.017 
0.009 
0.010 
0.009 
0.014 
0.017 
0.017 
0.013 
0.018 
0.010 
0.020 
0.020 
0.022 
0.010 
0.010 
0.010 
0.010 
0.009 
0.010 
0.009 
0.010 
0.010 
0.009 
0.021 
0.017 
0.021 
0.010 
0.010 
0.009 
0.020 
0.022 
0.021 
0.021 
0.010 
0.022 
0.011 
0.015 
0.011 
0.012 
0.027 
0.015 
0.037 
0.067 



0.881 
0.888 
0.944 
1.017 
1.061 
1.116 
1.230 
1.309 
1.353 
1.355 
1.348 
1.522 
1.511 
0.614 
0.668 
0.722 
0.768 
0.837 
0.875 
0.883 
0.860 
0.965 
1.100 
1.156 
1.121 
1.078 
1.234 
1.176 
1.225 
1.268 
1.510 
1.330 
1.342 
1.406 
1.441 
1.450 
1.513 
1.504 
1.549 
1.570 
1.393 
1.577 
1.625 



0.023 
0.013 
0.015 
0.013 
0.019 
0.025 
0.025 
0.019 
0.030 
0.016 
0.030 
0.035 
0.057 
0.014 
0.014 
0.014 
0.014 
0.013 
0.014 
0.013 
0.014 
0.014 
0.013 
0.030 
0.025 
0.028 
0.014 
0.015 
0.013 
0.029 
0.033 
0.031 
0.031 
0.017 
0.036 
0.021 
0.027 
0.025 
0.032 
0.064 
0.063 
0.045 
0.071 



0.956 
1.051 
1.039 
1.118 
1.262 
1.237 
1.424 
1.505 
1.550 
1.619 
1.701 
1.863 
1.999 
0.718 
0.786 
0.805 
0.844 
0.912 
0.959 
0.941 
0.994 
1.058 
1.233 
1.289 
1.257 
1.328 
1.373 
1.378 
1.634 
1.498 
1.901 
1.584 
1.600 
1.749 
1.815 
1.911 
1.934 
1.963 
2.137 
2.178 
2.285 
2.350 
2.626 



0.024 
0.013 
0.014 
0.013 
0.027 
0.024 
0.024 
0.018 
0.025 
0.014 
0.027 
0.027 
0.029 
0.015 
0.014 
0.015 
0.014 
0.013 
0.014 
0.013 
0.014 
0.014 
0.013 
0.029 
0.024 
0.029 
0.013 
0.014 
0.013 
0.028 
0.034 
0.030 
0.029 
0.013 
0.030 
0.015 
0.021 
0.015 
0.016 
0.035 
0.019 
0.045 
0.071 



Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 



Yes-sing 
Yes-sing 
Yes-sing 
Yes-sing 
Yes-bin 

Yes-sing 
Yes-sing 



Yes-sing 

Yes-sing 

Yes-sing 

Yes-sing 

Yes-sing 

Yes-sing 

Yes-sing 

Yes-bin 

Yes-sing 

Yes-bin 

Yes-sing 

Yes-bin 
Yes-sing 

Yes-bin 



Table 2. continued. 
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pi 01 6S 
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4-49 46 90 47^ 
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4-49 4Q ^ 840 
+*+Z *+V JJ.OHU 


IP 9 1 8 
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1 ft 949 
ID. Z^fZ 
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01 s 
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Mn 
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F3D437 


02 40 48.970 


+42 56 12.434 




16.022 


0.018 


0.647 


0.025 


o 


.748 


0.025 


No 




F3.0692 


02 40 57.994 


+42 49 44.614 




17.136 


0.010 


0.857 


0.016 


1 


.061 


0.015 


No 




F3^310 


02 40 44.085 


+42 53 58.059 




99.999 


9.999 


9.999 


9.999 


9 


.999 


9.999 






F4.0167 


02 4148.825 


+42 31 31.702 


JP223 


13.502 


0.021 


0.938 


0.026 


1 


.030 


0.031 


No 




F4.0404 


02 41 27.663 


+42 35 42.306 


JP131 


14.730 


0.010 


0.852 


0.014 





.932 


0.014 


Yes 


No 


F4.0515 


02 41 27.188 


+42 43 41.881 


JP 129 


15.111 


0.010 


0.903 


0.014 





.977 


0.014 


No 





o 

I 



a— Identifier format is XX_YYYY, where 'XX' are representations F3 and F4 for targets detected on M 34 western Field 3 or western Field 4 images, 
and 'YYYY' are internal catalog numbers. 

b— Astrometric data are derived from coordinate matches in the SuperCosmos survey catalogs (Hambly et al. 2001a, b, c). 
c- JP 96 identifiers are for photometric variables which are correlated with Jones & Prosser (1996) sources 
d- BVI photometric data are taken from the Indiana Group survey (c.f., §[3). 

e— Membership assignments are based upon whether a target star lies on (or close to) the V/B-V and/or V/V-I main sequence loci [see Fig., Q~) or whether a target star has a 
radial velocity consistent with the cluster's systemic velocity (Meibom et al. in prep; see also § 12. H and Meibom et al. 2006). 

/- B-V colours are determined using a field star V-I to B-V relationship (Caldwell et al. 1993), assuming Ey-i - 1.25xEb-v, and Eb-v-0-01 (Canterna et al. 1979). 
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Table 3. Photometric periods and ancillary data are presented, resulting from the light-curve analysis of variable stars detected in the Lowell M 34 
fields. 



Internal" 


Data 


Sigma ClippecF 


Period 


B-V„ 


Mem 


Identifier 


bpocns 


StdDev [mags] 


[days] 






T~"~i A 1 T O 1 in r \ f 

F3_0172-lIlVsV 


C ") /C 1 IC A 

53/54/59 


A AAOOO /A A 1 A 1 A /A A 1 1 AA 

0.00883/0.01010/0.01 190 


8.0 ± 1.0 


0.811 


Y;Ys 


r3_017o-lll Vslsv 


CO /C .1 /CO/f A 

53/54/58/59 


A A A T7C /A A/I /I C A /A A/11 1 1 /A A/C /I n 

0.04275/0.06460/0.04313/0.06427 


O "7 OA 1 A AAC 

2.780 ± 0.005 


A O 1 O 

0.818 


Y;Ys 


r3_02 15-111 Vslsv 


52/54/58/59 


A A 1 10 1 /A A1TWA A1 A OC\ /A AO AO A 

0.01 Izl/0.017zz/0.0143y/0.0z0z0 


HAH 1 A 1 1 

7.47 ± 0.1 1 


A OH A 

0.874 


Y;Ys 


b3_0258-WVsV 


CO /c^ /CO 

53/54/58 


A A1 IOO /A Al 1CA /A A1 A C7 

0.01 182/0.01250/0.01457 


1 1 A 1 1 A 

11.0 ± 1.0 


A A/17 

0.947 


Y;Ys 


F3_0306-nsIsV 


C O /C O /c a 

53/58/59 


A A 1 O O f /A A 1 O AA /A AO A O 1 

0.01886/0.01809/0.02081 


1 1 A O C A AAAT 

1.1435 ± 0.0007 


0.991 


Y;Yb 


T'O AT>A ITIT 7 ^7 

F3_0320-lIlVsV 


C O /C /I /C O 

52/54/58 


A AAO A O /A A 1 ATI /A A 1 A 1 A 

0.00898/0.01023/0.01430 


O C 1 1 A AA 1 

8.511 + 0.001 


1.046 


Y; 


F3_0383-1I1V 


52/54 


A AAT"7/^ /A AAA1H 

0.00776/0.00937 


9.5 + 0.4 


1.160 


Y;Ys 


ni A/| TA 1 \ 7 

r3_0430-lV 


C 1 

54 


A A 1 OCT 

0.01253 


1 A C 1 1 A 

10.5 ±1.0 


1 OO A 

1.239 


Y;Ys 


r 3 _04o9-lll V sis V 


52/54/57/57 


A A 1 O /l /I /A A1T/IC/A Al 7£0 /A AO O O C 

0.01344/0.01745/0.01768/0.02225 


1 O A 1 AO 

12.9 ± 0.3 


1 OOO 

1.283 


Y; 


1 ""l f\AOC 1T1\7 

r3_0485-lllV 


52/54 


A Al A/IA/A C\ 1 "> A "> 

0.01040/0.01393 


1 O A 1 1 A 

1Z.0 ±1.0 


1 O O C 

1.285 


Y; 


r3_0487-WV 


53/54 


A A 1 A7 /T /A A 1 IOO 

0.01076/0.01388 


OA 1 A H 

8.4 ± 0.6 


1 •"no 

1.278 


Y; 


01 r\tzr\r\ iti\/ 
r3_0o00-lllV 


CI /C /I 

53/54 


A A1 T/CI /A AOCAT 

0.01263/0.02503 


H A A 1 A AC 

/.44 + 0.05 


1 A CO 

1.452 


1 ; 


F3_0664-1I 


52 


A A 1 A^T 

0.01077 


10.5 ±1.0 


1.441 


Y; 


r4_013o-slsV 


60/60 


A A1 A CI /A A 1 ~7 A /I 

0.01457/0.01794 


O OH 1 A AO 

2.87 ± 0.03 


A C A A 

0.544 


Y;Ys 


C/1 A1 £A 1T1\7'„T„\7' 

r4_0 1 09-111 VslsV 


ccmo/^a/ca 
55/42/60/59 


A A1 17C/A A1 OAT /A A 1 A A A /A A1 /I OA 

0.01Z /5/0.01Z9 //0.010V4/0.014Z9 


*3 O 1 i A AO 

3.81 ± 0.08 


A CAO 

0.598 


Y; Ys 


F4_0194-slsv 


60/60 


A A A on /A A A O O A 

0.00817/0.00889 


O OA A AO 

2.20 ± 0.03 


0.652 


Y;Ys 


F4_0234-lDVsV 


54/55/60 


0.00692/0.00984/0.01014 


5.56 ± 0.07 


0.698 


Y;Ys 


T7/1 AO AO 1T1 T 7„T^\ 7 

r4_0303-lll VslsV 


jj/jj/oU/oU 


A AAOAA/A AAA/1 1 /A AAO T A /A AAATO 

0. OOoOO/0. 00941/0. 00o39/0.00V /o 


HA 1 A 

/.4 ± 0.3 


A HUH 

0. /of 


V. V.-. 

Y; Ys 


r4_03 17-111 Vslsv 


cc /cc /cn icc\ 

55/55/59/60 


A A 1 TA 1 /A AO 1 C A /A A 1 AO /A AO 1 ~7 

0.01701/0.02159/0.01698/0.02167 


H O 1 A 1 

6.8 ±0.1 


A O AC 

0.805 


Y;Ys 


F4_0327-lVsV 


55/60 


a a 1 /r\ a 1 /i a*t 

0.01471/0.01406 


7.6 ±1.1 


0.813 


Y;Ys 


r4_U33j-lll Vslsv 


CC /CC If- A //i A 

jj/jj/oU/oU 


A Al COO /A A1ACT/A A1 A 1 /A AO AO O 

0.0153o/0. 01953/0. 01o01/0.0z03z 


A OA 1 A AO 

0.89 ± 0.02 


A 7AA 
0. /90 


Y; Yb 


A A/1 CA 1T1\7~T~\7 

r4_0450-lll VslsV 


C /I /c c //; A Id A 

j4p5/oU/oU 


A AAA OC/A A1 1 01 /A A10/C/1 /A A1 I^O 

0.009o5/0.01 lol/0.01zo4/0.013oo 


O AA 1 A 1 

8.00 ± 0.21 


A OAC 

0.895 


Y; Ys 


T7/1 A/C/C-7 1T1\7„\7 

r4_Uoo7-lllVsV 


54/55/59 


A A1 IOO /A A1 OIA /A AO AO O 

0.01 Izz/0.01ol9/0.0z0zo 


1 .91 + 0.05 


1 AO A 

1.030 


Y;Yb 


T7/1 A/CAC IT 

r4_0695-ll 


C /I 

54 


A AAT /I O 

0.00748 


/I O 1 A C 

4.3 ± 0.5 


1 AO/C 
1.086 


Y;Ys 


r4_0 / 30-111 V sis V 


cc /cc /^A/cn 
55/55/00/59 


A AO A A A /A A/I/;1 A/A AO 1 AA /A A /I /I AC 

0.0Z944/0.040 19/0.03109/0.04095 


O A 1 i A 1 A 

0.41 + 0.10 


1 AC 1 

1.051 


y ; 


r?/i mii: iti\7„t„\7 

ir'4_0736-lll VslsV 


C C /C C ICC\ICf\ 

55/55/60/60 


A A 1 0"7 C /A AOT1 C/A AOAOO /A AO A 1 1 

0.01875/0.02715/0.02022/0.0291 1 


A HOU 1 A AA 1 

0.786 ± 0.001 


1 AAO 

1.008 


Y;Yb 


r~< a AOAO IT 7 T 

F4_0803-lvsI 


55/60 


A A 1 /11C /A Al 41 / 

0.01435/0.01416 


12.6 ± 0.7 


1.164 


Y;Ys 


T7A AOA/I 1T1\7„T„\7 

r4_0804-lll VslsV 


55/55/60/60 


A AO C/C O /A A A 1 CO /A AO AAO /A A/IOTA 

0.02568/0.04152/0.02448/0.04239 


A AC 1 A A A 

4.95 ± 0.04 


1.106 


1; 


f~> A aqCI 1T1T7 T T 7 

F4_0852-lIlVsIsV 


CC /CC /CA l/Zf\ 

55/55/59/60 


A A 1 C A 1- ! /A AO 1 ^ O /A A 1 ^ O ^ /A AO C tl A 

0.01547/0.02368/0.01624/0.02564 


A OHC\ 1 A AAO 

0.879 ± 0.002 


1.155 


Y;Yb 


r4_0942-WVsV 


C< /CC //C A 

54/55/60 


A AAA A O /A Al 1 AA /A AO T/TI 

0.00942/0.01300/0.02362 


A A 1 A O 

9.4 ± 0.3 


1 1 AO 

1.198 


Y; 


a 1 aa^ 1 ti \ 7 

F4_1007-1I1V 


55/55 


A AAT"7 /I /A A 1 1 m 

0.00774/0.01 177 


15 ± 2 


1.440 


Y; 


V A 1 A1A 1T1T7 T \ T 

F4_l 020-111 Vslsv 


CC /CC /ZTA/CA 

55/55/60/59 


A A 1 O A A /A A 1 O A H /A A 1 C O /A AO COO 

0.01344/0.01846/0.01568/0.02528 


8.6 ± 0.4 


1.260 


Y; 


r~< a 1 f\C C 1T1T7 T \ T 

F4_1055-lIlvsIsV 


C C IC C l£Lf \ If f\ 

55/55/60/60 


A A 1 ;in /r\ AO C y| A //"» r\oo 11 /A AO tl tl A 

0.01427/0.02542/0.02233/0.02664 


tl OA 1 A 1 C 

6.80 ± 0.15 


1.272 


Y; 


11/1-7 H 7 

r4_l 14 /-IV 


55 


A A 1 OOI 

0.01883 


A C O A 1 A AA 1 

0.589 ± 0.001 


1.336 


"V". 

1; 


T7A ITIf 1T1\/ 

r4_l jlj-lilV 


Cyl /C/l 
54/54 


A Al A /A AlOll 

0.0142y/0. 01913 


1 O O _i_ A C 

Iz.z ± U.j 


1 1"7 1 
1 . J / 1 


1; 


F4 1 357-111 V 


55/55 


01 197/0 01973 


9.4 ± 0.2 


1.380 


Y; 


PA 1 ^AQ 1T1V 




U.UlooO/U.UZj / 


11 ss + nn 

1 1 .0 J ± W. 1 J 


1 AAT, 
1 .44 J 


v- 




JJ/JJ 


U.UIt-^U/U.UZjOj 


ft s + s 

O.J X O.J 


1 4^4 

1 .4 J4 


v- 

1 > 


PA 1 Sfifi-ITIV 

I t.lJUU 111 V 


SS/S4 


01 Q70/0 0^?A8 


11 9 + 01 
1 1 .z in u. 1 


1 47Q 
1 .4 / 7 


Y- 


PA 1617 IV 
r i +_loi /-IV 


JZ 




/ . J ± O.J 


1 SOO 
1 . JOO 


Y- 

i) 


PA 1 7Q^-1T1V 

I / 7J 111 V 




07771/0 OAQfSfS 


7 C, + A f. 


1 ^9^ 

1 . JZ J 


Y- 


PA 1 09^ IV / 


AA 


A AAO 1 A 

u.uyziu 


1 9 7 _i_ 9 
1Z. / ± Z.O 


1 ^07 
1 .JO / 


I, 


r4_ZZZo-lF 


55 


A AO/CO 1 

O.OZool 


A A _l A C 

0.0 ± 0.5 


1 ccc 
1.555 


1; 


r?i aah 1 t \ 7 

F3_0071-slsv 


58/59 


A AOACO/A AOCAO 

0.02052/0.02502 


O A 1 A 1 AAA O 

3.014 ± 0.008 


0.977 


N; 


FTO AlTl 1\7„\7 

b3_0121-lVsV 


C A ICC\ 

54/59 


A AAOC1 /A AAO OI 

0.00851/0.00881 


OA 1 A O 

8.4 ± 0.Z 


1 1 "7 O 

1.178 


IN; 


b3_01o3-llslsV 


CO /CO /CA 

53/58/59 


A AAO ~\ A lf\ AAA /I O /A AAA 1 C 

0.008 14/0.00943/0.00935 


i/i 1 
14 ± 2 


A OO A 

0.830 


Y;N 


T70 A 1 £ C IT 


CT 

52 


A A 1 AO A 

0.01089 


1 1 1 
12 ± 1 


A HC\C 

0.795 


"\/, X T 

Y;JN 


O^RR ITIVcTcV 
r j_ojoo-in v sis v 




n m s^zt/n 091 A'\ /n nni s/n 09ziziR 

W.Wl JJ4/W.WZ141/W.W1 / 1J/W.WZ440 




QOS 
O. VOJ 


NT- 


F3.0413-lIlVsIsV 


53/54/58/59 


0.01640/0.02417/0.01909/0.02982 


5.45 ±0.08 


1.073 


Y;N 


F3.0437-lIlVsV 


53/54/59 


0.00859/0.01134/0.01563 


9.7 ± 0.4 


0.577 


N; 


F3.0692-1I1V 


53/54 


0.01886/0.02686 


0.3368 ± 0.0001 


0.787 


N; 


F3 .23 10-111 V 


52/54 


0.02667/0.02846 


2.991 ±0.008 






F4.01 67-111 Vsl 


54/55/59 


0.01321/0.01091/0.00961 


9.2 ± 0.9 


0.868 


N; 


F4.0404-lIlVsV 


55/55/60 


0.01099/0.01419/0.01461 


6.6 ±0.3 


0.782 


YN 


F4.0515-lIlVsV 


55/55/60 


0.00886/0.00957/0.01239 


8.2 ±0.1 


0.833 


N; 



a- Internal Identifiers are the same as those detailed in Table [2] Designations indicating exposure time and filter choice for a given 
light-curve are annotated. For instance, //represents long-I, for light-curves derived from long exposures in the I-filter. Thus, MVsIsV 
represents four independent light-curves created from image sequences of 400seconds, I-filter; 600seconds, V-filter; 40seconds, I- 
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filter & 60seconds, V-filter, respectively (namely, long-I, long-V, short-I & short- V). 

b— Number of photometric data points in each light-curve, derived from each target's image series, in a given filter and exposure 
time sequence. 

c— For each target's light-curve, in a filter+exposure combination, standard deviations of the differential magnitudes are provided. 
The time series files for each target were sigma-clipped to eliminate outliers. A A-magnitude for each light-curve can be determined 
by multiplying these standard deviations by 2 V2 (strictly speaking this is only true for perfectly sinusoidal variations). 
d- E B _y=0.07 is adopted; Canterna et al. (1979) 

e— Abridged cluster membership criteria are reproduced from Tabled with format Photometric; Kinematic (s-single, b-binary). 
/- B-V colours are determined using a field star V-I to B-V relationship (Caldwell et al. 1993), assuming Ey-i = 1.25xEb-v, and 
E B _ V =0.07. 
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Table 4. A comparison of photometric periods for three M 34 
stars common to the Lowell campaign, the Irwin et al. study and 
the Meibom et al. (in prep) study (c.f., Fig.,|2]and § I4.4| |. 



ID 


Lowell 


Irwin et al. 


Meibom et al. 




[days] 


[days] 


[days] 


F3.0413 


5.45 ± 0.08 


2.458 


5.47 


F4.0234 


5.56 ± 0.07 


0.822 


5.47 


F4.1147 


0.589 ± 0.001 


1.429 


0.592 
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Table 5. Gyrochronological ages for I-sequence stars in M 34 



Identifier 


JP" 


#of 


(B-V)o 


Measured 


Gyro Age r 


Fitted d 


Period^ 




no. 


periods 




Period 




Period 


Variance 






derived 6 




[days] 


[Myr] 


[days] 


[days 2 ] 


F4.0136 


JP213 


2 


0.544 


2.87 ± 0.03 


178.5 ± 3.5 


2.99 


0.0144 


F4.0169 


JP50 


4 


0.598 


3.81 ±0.08 


170.0 ±6.7 


4.08 


0.0729 


F4.0234 


JP 148 


3 


0.698 


5.56 ± 0.07 


188.4 ±4.4 


5.63 


0.0049 


F4.0303 


JP227 


4 


0.767 


7.4 ±0.3 


244.1 ± 18.6 


6.5 


0.81 


F4.0317 


JP224 


4 


0.805 


6.8 ±0.1 


183.8 ±5.0 


7.0 


0.04 


F4.0327 


JP 199 


2 


0.813 


7.6 ± 1.1 


220.9 ± 59.8 


7.1 


0.25 


F3.0172 


JP49 


3 


0.811 


8.0 ± 1.0 


244.6 ± 57.2 


7.0 


1.00 


F3.0215 


JP41 


4 


0.874 


7.47 ±0.11 


180.4 ±5.0 


7.74 


0.0729 


F4.0450 


JP 113 


4 


0.895 


8.00 ± 0.21 


194.5 ± 9.6 


7.96 


0.0016 


F3.0258 


JP289 


3 


0.947 


11.0 ± 1.0 


313.1 ±53.3 


8.5 


6.25 


F4.0667 


JP52 


3 


1.030 


7.97 ± 0.05 


145.0 ± 1.7 


9.28 


1.7161 


F3.0320 


JP 172 


3 


1.046 


8.511 ±0.001 


159.3 ±0.1 


9.429 


0.84272 


F4.0730 


JP 18 


4 


1.051 


8.41 ±0.10 


154.4 ±3.4 


9.47 


1.1236 


F3.0383 




2 


1.160 


9.5 ± 0.4 


162.2 ± 12.8 


10.4 


0.81 


F4.0803 


JP229 


2 


1.164 


12.6 ± 0.7 


273.5 ± 28.5 


10.5 


4.41 


F4.0942 




3 


1.198 


9.4 ±0.3 


150.4 ± 9.0 


10.7 


1.69 


F3.0430 




1 


1.239 


10.5 ± 1.0 


174.8 ±31.2 


11.1 


0.36 


F3.0485 




2 


1.285 


12.0 ± 1.0 


211.3 ±33.0 


11.4 


0.36 


F3.0469 




4 


1.283 


12.9 ±0.3 


242.5 ± 10.6 


11.4 


2.25 


F4.1315 




2 


1.371 


12.2 ±0.5 


196.4 ± 15.1 


12.1 


0.01 


F4.1357 




2 


1.380 


9.4 ± 0.2 


119.3 ±4.8 


12.2 


7.84 


F3.0664 




1 


1.441 


10.5 ± 1.0 


137.2 ±24.5 


12.6 


4.41 


F4.1369 




2 


1.443 


11.85 ±0.13 


171.7 ±3.5 


12.61 


0.5776 


F4.1566 




2 


1.479 


11.2 ± 0.1 


148.8 ± 2.5 


12.9 


2.89 


F4.1007 




2 


1.440 


15.0 ±2.0 


267.8 ± 66.8 


12.6 


5.76 


F4.1925 




1 


1.507 


12.7 ±2.0 


183.0 ±53.9 


13.1 


0.16 



a— JP identifier from Jones & Prosser (1996). 

b- Number of filter/exposure combinations yielding photometric period derivations for a given star. 

c- Gyro ages determined using Meibom et al. (2009) I-sequence colour function coefficients, viz a = 0.770 ± 0.014, b = 0.553 ± 0.052 & 
c = 0.472 ± 0.027, and an age function exponent n = 0.5344 ± 0.0015 (see § |5~TT >. 

d— Assuming a mean gyrochronology age for M34 of 193Myr (the mean of the ages calculated in column 6 - see S I5.lt and the Meibom et al. 
(2009) coefficients, with n = 0.5344 ± 0.0015, we use intrinsic (B-V) colours and Equ.,1 to derive a fitted period for each M 34 I-sequence star. 
e— Each star's contribution to the period variance [(Measured Period - Fitted Period) 2 ] is noted, where the total variance is 43.73 days 2 , with a 
mean of 1.68 days 2 . 
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B-V V-l 

Fig. 1. Colour magnitude diagrams, in B-V and V-I colours, are presented for the entire photometry catalogue of the M34 BVIc 
survey conducted by the Indiana Group. Solid lines represent D' Antona & Mazzitelli (1997) theoretical isochrones, for a cluster age 
of 250 Myr, a distance of 470pc, and assuming Ey-i - 1 .25 x Eb-v, and Eb-v-0-01 (Canterna et al. 1979). Red squares, highlighting 
photometric and radial velocity members of the cluster (solid symbols), and photometric cluster membership only (open symbols) 
are plotted for M 34 stars for which we have derived photometric periods during the Lowell campaign. Blue circles depict those 
other periodic photometric variables in the field of M 34 which are photometric and/or radial velocity non-members of the cluster. 
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Period [days — Lowell Campaign] 



Fig. 2. Photometric periods for M 34 stars common to our Lowell campaign and the Irwin et al. study are shown. Periods determined 
from the Lowell campaign, using photometric light-curves in up to four (4) different filter/exposure time combinations, are depicted 
by the filled-square, circle, triangle and asterisk symbols (see Tables [2] & 0. We include error bars on our Lowell periods, which 
for many targets are smaller than the symbol sizes. No error bars are available for the Irwin et al. periods. The solid line represents 
equality between the two studies, and is not a fit to the data. The dotted line represents phase-halving while the dashed lines represent 
several period-aliasing loci, respectively. About one-third of the sample have equal period determinations (within the errors). Of the 
remainder, all but two of the Irwin et al. periods are shorter than ours. See text for possible reasons for the discrepancy (see § !4.4t . 
The three flagged stars [with boxed symbols] are discrepant ones where our periods have been confirmed by a forthcoming Meibom 
et al. (in prep) study. 
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Fig. 3. Photometric period distributions, as a function of intrinsic colours (B-V) c and (V-I) G , are plotted for solar-type stars in NGC 
1039 (M34). We adopt Eb_v=0.07 (Canterna et al. 1979), and Ey-i - 1.25x Eb-v- Period symbols are the same as those depicted 
in Fig., |2] Error bars are included, and are for some targets smaller than the symbol sizes. Approximate spectral types for a given 
colour are annotated on the upper abscissae (from Zombeck 1990). 
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(B-V)o (B-V)o 



Fig. 4. In the left-hand panel, the rotation period distribution (in days) of the ^ 1 35 Myr M 35 cluster is plotted against intrinsic (B- V) 
colour (data from Meibom et al. 2009); stars which are both photometric and kinematic members of the cluster are plotted with larger 
symbols than stars having photometric membership alone. Over-plotted are 135 Myr interface [solid line] and convection [dashed 
line] gyrochronology loci from the calibration by Meibom et al. (2009). In the right-hand panel, similar colour-period data are 
plotted for the 600 Myr Hyades cluster (data from Radick et al. 1987; Prosser et al. 1995b). Under the paradigm of gyrochronology, 
specific regions of the colour-period distributions are identified with C [convective], g [gap] and I [interface] symbols. 
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Fig. 5. The colour-period distribution for all photometric and/or kinematic members of M 34 (see Table [3} is plotted in the left 
hand panel. Over-plotted are gyrochronology Interface and Convection sequence loci (solid & dashed lines respectively) for the 
135 Myr M35 cluster (Barnes 2003, Meibom et al. 2009). Regions of the colour-period diagram where approximately 200 Myr 
stars, classified as Interface, gap and Convective status in gyrochronology space, are identified with I, g or C symbols respectively. 
Candidate Interface sequence stars in the M 34 cluster, lying above the Meibom et al. M 35 Interface sequence, are identified as 
I-symbols in the right-hand panel; see §[5]for details. 
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Fig. 6. This figure, very similar to Fig., [5] shows the colour-period data for our M34 candidate I-sequence stars. The additional 
solid line in the right-hand panel represents the mean gyrochronology age of the I-stars in M 34 (193 Myr - see Table[5]l. This new 
gyrochronology I-sequence is calculated using updated colour-function coefficients (see Equ.,1-3 and Meibom et al. 2009) together 
with our new derivation of the g(t) exponent (n = 0.5344 - see Equ., 3). 
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Appendix A: Cluster Membership of M 34 Using 
Extant Kinematic Measurements 

There exist two extant kinematic studies of the M34 cluster 
which can assist us further in establishing cluster membership 
for our gyrochronology sample. The first, by Jones & Prosser 
(1996), relies upon 2-d space motions in the form of proper mo- 
tions, although their high-fidelity membership probabilities are 
magnitude limited to about 14.5. While JP96 proper motion 
results do extend to fainter magnitudes, their reliability reduces. 
That is to say, a fainter M 34 star with a 00% proper motion 
probability may indeed still be a bona fide member of the clus- 
ter, insomuch as Poisson-limited proper motion errors could lead 
to erroneously lower proper motion probabilities. 

The second, by Jones et al. (1997), employing 1-d heliocen- 
tric radial velocities, can be used to establish cluster member- 
ship, or in the very least, establish cluster non-membership (or 
binarity). Unfortunately neither membership probabilities nor 
velocity errors were reported by Jones et al., however we can 
exploit their measurements in our membership assessment (see 
below). Interestingly, in the same manuscript, Jones et al. report 
Li I 6708A equivalent widths [EWs] for their target stars. Even 
though the detection of lithium in solar type stars is not a req- 
uisite for cluster membership, its absence is almost ubiquitous 
among the generally older, Galactic field stars. The very pres- 
ence of a substantial lithium line in solar-type stellar spectra is 
in itself indicative of youth (<< 1000 Myr; e.g., Soderblom et 
al. 1993b, James & Jeffries 1997, James et al. 2006). 

For each of the M34 stars in our sample, we have col- 
lated these additional membership criteria published by JP96 and 
Jones et al. (1997) in Table IA.fi Its Lowell identifier is listed in 
column 1 in concert with its JP96 identifier in column 2. Each 
star's photometric and radial velocity membership flag, repro- 
duced from Tables [2] & [3] are cited in columns 3 & 4 respec- 
tively, whereas a flag confirming the detection of lithium 6708 A 
in its spectrum (from Jones et al. 1997) is listed in column 5. 
Corresponding Jones et al. radial velocities are detailed in col- 
umn 6, as well as their associated membership probabilities in 
column 7 (see below and Fig. IA. 11 1. V-band magnitudes of each 
star are listed in column 8 and finally, cluster membership prob- 
abilities based on the proper motions reported in JP96 are repro- 
duced in column 9. 

In the first instance, we note that all stars in common between 
our gyrochronology sample and the Jones et al. sample have a 
significant lithium detection, indicative of stellar youth. These 
stars are therefore unlikely to be field star interlopers, which 
when combined with the photometric membership criterion, can 
safely be considered as bona fide cluster members. In the second, 
the radial velocity data from Jones et al. study do seem to cluster 
around -10km s -1 , however, specific membership assignments 
are difficult to assess due to the associated scatter in the indi- 
vidual measurements. In order to understand whether the scatter 
in these radial velocities is due to a true dispersion in the data 
or is an admixture of single and binary member measurements, 
together with cluster non-members, we have constructed a his- 
togram for all radial velocities published in the Jones et al. sur- 
vey. This histogram, plotted in Fig. IA.ll clearly shows a peak in 
the distribution representing the cluster's systemic heliocentric 
radial velocity. In order to investigate the dispersion about the 
central peak, we fit an unweighted Gaussian function to the ve- 
locity histogram, noting that the fit does not include background 
contamination due to binarity or cluster non-members, and as- 
sumes that the sample is complete. The Gaussian function is cen- 
tred on -12.49 km s _I with a lcr width of 2.32 km s _1 , which 



we can exploit to determine a membership probability (reported 
in column 7 of Table IA.U for each individual star based on its 
radial velocity and these Gaussian fit parameters. 

A.1. Results: 

Analysis of the additional membership data reported in 
Table IA.fl yields few surprises. Of all the Lowell periodic vari- 
ables identified as photometric and/or RV members in Table [2] 
only two of them, F3_0172 & F3_0215, have a zero probability 
of being proper motion members according to the JP96 survey. 
Interestingly however, their V-magnitudes lie right at the point 
at which the JP96 proper motions begin to become less reliable 
due to Poisson-error limits on their photographic plate measure- 
ments, which brings their validity into question. Furthermore, 
we have noted in § !5.2l that the gyrochronology sample may in- 
deed suffer from contamination due to one or two cluster non- 
members masquerading as bone fide members, but we argue that 
their contribution to the period variance, and hence gyro errors, 
is small. In fact, F3_0172 & F3_0215 contribute a combined pe- 
riod variance of 0.79 days 2 , which is only a 1.8% effect (see 
Table|5]l. In spite of their low contribution to the period variance, 
we remain in the uncomfortable position of choosing whether to 
include these photometric and kinematic members of the cluster 
in our gyrochronology sample, or to exclude them on the ba- 
sis of their proper motion membership probabilities. In order to 
further assess their cluster membership status, we have recently 
obtained high resolution optical spectra of these two stars using 
the Canada France Hawai'i Telescope [CFHT], the analysis of 
which we discuss below (see § IA.2b . 

Finally, all bar one of the Jones et al. (1997) stars are radial 
velocity members of the cluster, which correlates well with our 
own radial velocity data. The one Jones et al. star, F4_0136, that 
is formally a cluster non-member based on a Gaussian probabil- 
ity fit to their radial velocity data plotted in Figure lATl actually 
shows up as a single cluster star member in the long term syn- 
optic velocity survey of Meibom et al. (in prep - see also § 12. U . 
Curiously, over the 9-epochs of observation covering 2.5-years, 
that Meibom et al. have for this star, its varies about their average 
radial velocity of -8.0 km s _1 by only 0.68 km s _1 . Interestingly, 
a Gaussian fit to all radial velocities for the 70 single and binary 
M 34 stars for which Meibom et al. have obtained results yields 
a cluster systemic velocity of -7.59 ± 1.02 km s _1 . Assuming 
the Jones et al. velocity is not in error, this star may be a binary 
member of the cluster, albeit with either a long-period orbit or a 
considerably eccentric one. 

A.2. Cluster Membership Status for stars F3.01 72 & 
F3.0215: 

Two stars exhibiting variability in our differential photometric 
survey of the M 34 cluster, for which we derive periods, namely 
F3_0172 & F3_0215, present us with somewhat of a conun- 
drum. While these stars have photometric and radial velocity 
properties consistent with cluster membership, with relatively 
short photometric periods indicative of youth (compared to the 
Galactic field), and result in gyrochronology ages appropriate 
for a 200 Myr group of stars, they possess proper motion vec- 
tors incongruent with the remainder of the cluster. In an attempt 
to establish or refute genuine cluster membership for these two 
objects, we have recently observed them at high resolution us- 
ing the fibre-fed, bench-mounted ESPaDOns echelle spectro- 
graph, located in a Coude-like instrument chamber in the CFHT 
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Observatory. The primary goal of these observations is to de- 
tect lithium at 6708A, and measure its equivalent width in these 
stars, thereby confirming their relative youth and increasing their 
probability of being bona fide cluster members. 

During the evening of 27 January 2010, high resolution 
(R^ 60, 000) spectra were acquired for stars F3.0172 & F3.0215 
using the fibre-fed ESPaDOns spectrograph, in service with 
the 3.6-m Canada France Hawai'i Telescope located on top 
of Mauna Kea, Hawai'i, USA. The ESPaDOns spectrograph 
consists of a 79 lines mirT 1 echelle grating imaged onto a 
2048 x 4608 EEV CCD detector, having 13.5 /jm square pix- 
els, with photon input delivered by separate 100 fim (1.6 arcsec) 
diameter sky and target fibres. This set-up yields a FWHM of 
cross correlated ThAr arc lines of 0.271A at 6700A, and a com- 
plete wavelength range of 3699 — > 10481 A. Using this set-up, 
F3_0172 & F3_0215 were observed for exposures totalling 2000 
& 2140 seconds respectively, resulting in spectra with S/N^ 20 
at 6700A. 

For each of the targets, we exploit our CFHT ESPaDOns 
spectra in order to measure their heliocentric radial velocity and 
their Li I 6708 A EW. The removal of CCD instrumental effects, 
as well as the extraction of wavelength-calibrated spectra, have 
been achieved by two independent data reduction methodolo- 
gies. In the first, we performed the bias-subtraction, flat-fielding, 
spectral order tracing, optimal extraction and wavelength cali- 
bration using standard IRAF procedures. In the second, we used 
the direct output from Libre-ESpRIT (Donati et al. 1997), the 
dedicated pipeline software for the ESPaDOns spectrograph. A 
cross-match of radial velocities and EWs for each target from 
the two reduction methods reassuringly yields consistent results 
to within ±0.1 km s and ±3mA respectively. 

Heliocentric radial velocities were derived, relative to the 
well-exposed standard star HD 32963, in the spectral region of 
the Mg triplet lines (5104 -> 5176A). Cross correlation of HD 
32963 with other radial velocity standard stars observed during 
the CFHT programme shows that the zero-point in placing our 
velocities onto the standard system is ^0.2 km s _1 . The relatively 
low S/N of our target spectra results in radial velocity precision 
errors of ^0.5 km s _1 . 

In order to measure equivalent widths for the Li 1 6708 A line 
in our targets, each spectrum between 6586 — > 671 1 A, was nor- 
malized using continuum fitting after spectra extraction. Each 
EW was calculated using both the direct integration and the 
Gaussian fitting methods, whose values were within a few per- 
cent of each other. In the lithium region, Li I EWs include con- 
tributions from the small Fe i+CN line at 6707.44 A, leading to 
measured Li I EWs which are representative of a slightly (10-20 
mA) over-estimated photospheric Li presence. Soderblom et al. 
(1993b) report that this Fe line blend has an EW = [20(B-V) 
- 3]mA, determined through an empirical relationship for main 
sequence, solar-type stars. For each target star, we removed the 
Fe I line contribution before transforming Li I EWs into abun- 
dances, N(Li) - on a scale where log N(H)=12, using the ef- 
fective temperature-colour (B-V) relation from Soderblom et al 
(1993c), and the curves of growth presented in Soderblom et 
al. (1993b). Data products, radial velocity and lithium measure- 
ments, from the CFHT spectroscopic observations of F3_0172 & 
F3JJ215 are presented in Table lA.2l 

CFHT radial velocities of F3J3172 & F3J3215 are consis- 
tent with cluster membership of M 34 irrespective of whether 
we compare their individual values to the Jones et al. (1997) 
sample or to the Meibom et al. (in prep) one. For the Jones et 
al. sample, their kinematic membership probabilities are non- 



zero although they are quite low at 08% & 03% (for F3.0172 
& F3_0215 respectively). Their membership probabilities are far 
more convincing when compared to the Meibom et al. sample 
(with respective values of 79% & >99%). 

The lithium content of both F3_0172 & F3_0215 is substan- 
tial indicating that these stars are relatively young compared to 
the general Galactic field, whose solar-type stellar content is typ- 
ically old and has had sufficient time to have proton-burned con- 
siderable fractions of its natal lithium. In order for these stars to 
be judged as likely members of the M 34 cluster, not only must 
they contain lithium in their atmospheres, but it must quantita- 
tively fit into the mass-dependent lithium abundance distribution 
for the cluster. In Figure lATZl we plot the lithium abundances for 
F3_0172 & F3J3215 that we have determined from our CFHT 
spectra in concert with the extant mass-dependent lithium dis- 
tribution for M 34 stars (data taken from Jones et al. 1997, who 
employed identical (B-V) -Temperature, Fe I line correction and 
lithium curves of growth as we have used). It is clear that both 
F3_0172 & F3_0215 have lithium abundances consistent with 
the mass-dependent lithium distribution of the M 34 cluster, and 
must be considered lithium abundance members. 

JP96 proper motion vectors for F3_0172 & F3JJ215 indi- 
cate cluster non-membership. However their V-magnitudes ren- 
der JP96 2-d kinematic membership probabilities questionable, 
because at or about this magnitude, JP96 proper motion accu- 
racy and precisions begin to have strong dependencies on the 
Poissonian errors of their centroiding measurements. In consid- 
eration that these two stars are both photometric and kinematic 
members of the cluster, have photometric periods consistent with 
the remainder of the cluster's distribution, and have measured 
lithium abundances which lie right along the trend of the mass- 
dependent lithium distribution for the cluster, they are probable 
M 34 cluster members and we retain them in our gyrochronology 
sample. 
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Table A.l. Cluster membership assessments for Jones & Prosser (1996) stars in our period sample. 



Internal" 


JP96* 


Phot" 


RV" 


Li Menf 


RV J97 r 


RV Mem" 


V" 


PM Prob c 


Identifier 




Mem 


Mem 


J97 


[km/s] 


J97 


[mag] 


JP96 


F3_0172 


JP49 


Y 


Ys 








14.584 


00% 


F3_0176 


JP 167 


Y 


Ys 






... 


14.659 


44% 


F3_0215 


JP41 


Y 


Ys 








14.918 


00% 


F3.0258 


JP289 


Y 


Ys 


Y 


-10.8 


Y (47%) 


15.236 


76% 


F3.0306 


JP265 


Y 


Yb 


Y 


-17 


Y ( 5%) 


15.503 


53% 


F3.0320 


JP 172 


Y 




Y 


-9.8 


Y (25%) 


15.664 


03% 


F4.0136 


JP213 


Y 


Ys 


Y 


-5.7 


N ( 0%) 


13.151 


94% 


F4.0169 


JP50 


Y 


Ys 








13.510 


04% 


F4.0194 


JP 133 


Y 


Ys 


Y 


-11.8 


Y (77%) 


13.662 


96% 


F4.0234 


JP 148 


Y 


Ys 








13.959 


05% 


F4.0303 


JP227 


Y 


Ys 


Y 


-15.4 


Y(21%) 


14.366 


92% 


F4.0317 


JP224 


Y 


Ys 


Y 


-9.4 


Y (18%) 


14.410 


91% 


F4.0327 


JP 199 


Y 


Ys 


Y 


-9.0 


Y(13%) 


14.483 


85% 


F4.0335 


JP 158 


Y 


Yb 


Y 


-9.9 


Y (26%) 


14.546 


86% 


F4.0450 


JP 113 


Y 


Ys 


Y 


-13.5 


Y (67%) 


14.960 


84% 


F4.0667 


JP52 


Y 


Yb 








15.623 


26% 


F4.0695 


JP 197 


Y 


Ys 








15.643 


05% 


F4.0730 


JP 18 


Y 










15.741 


36% 


F4.0803 


JP229 


Y 


Ys 


Y 


-9.9 


Y (26%) 


16.025 


46% 


F3.0121 


JP212 


N 










13.741 


00% 


F3.0163 


JP258 


Y 


N 








14.440 


00% 


F3.0165 


JP244 


Y 


N 








14.419 


00% 


F3.0388 


JP 174 


N 










15.952 


00% 


F3_0413 


JP218 


Y 


N 








16.242 


17% 


F4.0167 


JP223 


N 










13.502 


00% 


F4.0404 


JP 131 


Y 


N 








14.730 


00% 


F4.0515 


JP 129 


N 










15.111 


00% 



a- Data taken from Tablesf2]& [D 

b— JP identifier from Jones & Prosser (1996). 

c— Cluster membership based on radial velocities and lithium 6708A equivalent widths, taken from Jones et al. (1997) [J97]. 

d- Radial velocity membership assessments are based on a Gaussian fit to the entire Jones et al. (1997) radial velocity sample (see also Fig. lA. It . 

e- Proper motion membership probabilities are taken from JP96. 
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Table A.2. CFHT spectroscopic data products for M 34 candidate members F3_0172 & F3_0215. 



Internal" 


JP96* 


(B-V)° 


Temp c 


HJD rf 


RV 


Lil + Fel e 


Fel' 


Li I 6708 A* 


N(Li) A 


Identifier 






[K] 


[(+2450000)] 


[km s~'] 


EW [mA] 


EW [mA] 


EW [mA] 




F3_0172 


JP49 


0.811 


5139 


5223.713 


-8.4 


115 


13 


102 


2.003 


F3_0215 


JP41 


0.874 


4946 


5223.712 


-7.4 


153 


14 


139 


1.985 



a- Data taken from Tables[2]& E] 

b— JP identifier from Jones & Prosser (1996). 

c— Temperature determined from (B-V)„ relationship in Soberblom et al. (1993c - equ. 3) 
d- HJD - heliocentric julian date 

e— Equivalent width measurement, by Gaussian fitting, of the blended Li I and Fe I 6708A lines. 

/- Fel 6707.441 A line contribution, determined by an empirical (B-V) colour fit (see Soderblom et al. 1993b). 

g- Li I 6708A minus the Fe I 6707.441A contribution. 

h— Lithium equivalent widths are transformed into abundance [N(Li)] - on a scale where log N(H)=12, using a curve of growth in Soderblom et 
al. (1993b). 



David J. James et al.: The Gyrochronology Age of M 34, Online Material p 5 



c 

"_Q 



<J) 

E 



o 



lO 



o 




-20 



5 







■5 







Heliocentric RV [km s ] 



Fig. A.l. Heliocentric radial velocities, in histogram form using 3 km s _1 bins, are plotted for M 34 stars using velocity data detailed 
in Jones et al. (1997). The red solid line depicts an unweighted Gaussian fit to the data, with a Gaussian centre of -12.49 km s _1 and 
a sigma of 2.32 km s~' . 
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Fig.A.2. Logarithmic lithium abundances N(Li) = 12 + log(Li/H) versus effective temperature are plotted for stars in the M34 
cluster (data taken from Jones et al. 1997). Red filled squares represent abundances for F3_0172 & F3_0215 (JP 49 & JP 41), 
determined using data reported in Table I A. 2l with error parallelograms based on temperature errors of 100 K and 15% equivalent 
width uncertainties. 
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Appendix B: A Period vs Equatorial Velocity 
Analysis 

A comparison of photometric period and projected equatorial ro- 
tation rate for M 34 stars is plotted in Fig., IB. II For the most 
part, photometric periods are determined from the Lowell cam- 
paign described in this manuscript, except for a few cases where 
Irwin et al. periods are employed where Lowell ones do not ex- 
ist. Spectroscopic velocities for stars in common to the Lowell 
campaign and Irwin et al. survey are obtained from Jones et 
al. (1997). Three loci are also shown in the figure, represent- 
ing equality between period and v sin i data for GO, KO & MO 
dwarfs, where sin i = n/4 is assumed for each. 

There are three notable characteristics to the plot. First, the 
clustering of data-points for stars with periods 6.5 — > 9 days 
and v sin / ^ 10- 14 km s _1 lie to the right of the GO dwarf locus. 
Given that the M 34 stars having photometric periods are late-F 
to early M-dwarfs, these data-points are most likely indicative of 
those stars whose sin / values are > tt/4. Second, there are four 
stars with periods < 3 days and 10 <vsin/ < 25 km s _1 , whose 
period, v sin i data place them considerably below even the MO 
dwarf locus in the diagram. Assuming that these stars are bona 
fide cluster members, whose period and vsin/ values are cor- 
rectly determined, these stars must be high inclination systems 
(sin i < 7i / 4; i.e., becoming more pole-on), whose true equatorial 
velocities are considerably higher. 

Finally, there is one M 34 star (F3_0258) whose period de- 
terminations, from our Lowell data and by the Irwin et al. study, 
are seriously discrepant (see also Fig.,|2j. The periodicity for this 
star was detected in three filter/exposure time observations dur- 
ing our Lowell campaign. It is interesting to note that the Irwin 
et al. period is almost half that of our Lowell one, and their lower 
value could be due to phase aliasing, power leakage in the power 
spectrum as a consequence of their shorter observing window or 
multiple spot groups on the surface on the star during its obser- 
vation. In any case, either both period determinations are incor- 
rect, or either value from the Lowell campaign or the Irwin et al. 
study is in error. 

If we assume that one of the periods for this star is correct, 
we can make some predictive statements as to its equatorial ve- 
locity. With an intrinsic B-V colour of 0.95, its spectral type 
on the main sequence would be K2 or K3, placing it close to 
the central locus of the three plotted in Fig.. IB. 1 1 (assuming the 
star is inclined 45° to the line-of-sight). This scenario is consis- 
tent with its Irwin et al. period of 6.655 days. Conversely, if the 
Lowell period of 11.0 days is correct for this object, and it is 
a single member of the cluster lying on the main sequence, its 
inclination angle must be higher than 45° (sin/ > tt/4), whose 
appearance is more face-on to the line-of-sight. Hopefully, more 
extensive photometric monitoring of this star will reveal its true 
nature. 
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Fig. B.l. Projected equatorial velocities (v sin i) are plotted for those M 34 stars having corresponding photometric periods derived 
during the course of the Lowell campaign. V sin ;' data are taken from Jones et al. (1997). For those cases where there is no Lowell 
campaign period available, similar data from the Irwin et al. (2006) sample are presented. One exception is shown for the star 
F3_0258, linked by a straight solid line, whose Lowell campaign and Irwin et al. periods are seriously discrepant (11.0 ± 1.0 and 
6.655 days respectively). Three loci are depicted representing equal period-vs-v sin i relationships (assuming sin; = n/A) for GO, K0 
and MO dwarfs (solid, dashed and dotted lines respectively). In order of descending mass, stellar radii of R/R Q = 1-08, 0.81 & 0.61 
were employed (Gray 1992). 
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Appendix C: Lowell Light-curves for 

Photometrically Variable Stars in the Field of 
M34 
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Fig. C.l. Period-phased photometric lightcurves for M 34 variables, derived from short-V observations. 
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Fig. C.2. Period-phased photometric lightcurves for M 34 variables, derived from short-I observations. 
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Fig. C.3. Period-phased photometric lightcurves for M 34 variables, derived from long-V observations. 
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Fig. C.4. Period-phased photometric lightcurves for M 34 variables, derived from long-I observations. 
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